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Abstract In diffusive Cu and Au quantum wires at finite transport voltage U the 
non-equilibrium distribution function f(E, U) exhibits scaling behavior, 
f(E,U) = f(E/eU), indicating anomalous energy relaxation processes 
in these wires. We show that in nonequilibrium the Kondo effect, gen- 
erated either by magnetic impurities (single-channel Kondo effect) or 
possibly by non-magnetic, degenerate two-level systems (two-channel 
Kondo effect), produces such scaling behavior as a consequence of a 
Korringa-like (pseudo)spin relaxation rate oc U and of damped pow- 
erlaw behavior of the impurity spectral density as a remnant of the 
Kondo strong coupling regime at low temperature but high bias. The 
theoretical, scaled distribution functions coincide quantitatively with 
the experimental results, the impurity concentration being the only ad- 
justable parameter. This provides strong evidence for the presence of 
Kondo defects, either single- or two-channel, in the experimental sys- 
tems. The relevance of these results for the problem of dephasing in 
mesoscopic wires is discussed briefly. 



1. INTRODUCTION 

Recently, it has become possible to determine experimentally [1, 2] 
the distribution function of quasiparticles, f x (E,U), in dependence of 
the particle energy E in a metallic, diffusive nanowire in stationary non- 
equilibrium with a finite transport voltage U applied between the ends 
of the wire. The measurements were done by attaching an additional 
superconducting tunneling electrode at a position x along the wire, so 
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that f x (E, U) could be extracted from the tunneling current 



jt = h No I dE [f*( E ' U ) ~ f°( E + eU t )]N BC s(E + eU t ) . (1.1) 



Here Nbcs is the peaked BCS density of states in the superconducting 
electrode, N Q and 7 are the approximately constant density of states in 
the wire and the tunneling matrix element, respectively, and Ut denotes 
the voltage between the wire and the tunneling electrode. f°{E) = 
l^ e E/k B T _|_ -g - g Fermi distribution at temperature T. 

As expected for a diffusive wire with elastic impurity scattering [3], 
f x (E,U) has a double-step shape, corresponding to the two chemical 
potentials at E = and E = eU at the ends of the wire. In addition, 
the steps are rounded, indicating inelastic scattering processes in the 
wire. In Cu [1] and Au [2] wires, the observed rounding is such that the 
quasiparticle distribution obeys a remarkable scaling property when U 
exceeds a certain low-energy scale, eU ~ eU ~ O.lmeV: f x (E,U) = 
f x {E/eU) [1]. One may gain phenomenological insight into the nature 
of the inelastic relaxation processes by observing [1] that the scaling 
property of f x (E,U) implies that the equation of motion of f x (E,U), 
the quantum Boltzmann equation, and consequently the inelastic single- 
particle collision rate 1/r are scale invariant as well. Assuming a (yet to 
be determined) two-particle potential V(e) with energy transfer e, 1/r 
is in 2nd order perturbation theory given as 



where F is a combination of distribution functions f x ensuring that 
there is only scattering from an occupied into an unoccupied state. The 
experimental finding that f x and hence 1/r depend only on the dimen- 
sionless energies, as indicated in Eq. (1.2), implies that in the term on 
the right-hand side only dimensionless energies (normalized to eU) oc- 
cur, i.e. V(e) oc 1/e in the scaling regime \E\ ~ k B U Q ~ eU. It follows 
that within 2nd order perturbation theory (Eq. (1.2)) the energy relax- 
ation rate has a logarithmic energy dependence, 1/t(E) oc \n(E/eU ), 
non-vanishing at the Fermi energy. This anomalous behavior has raised 
considerable interest in the energy relaxation measurements [1], espe- 
cially because of the possible relation of the apparently non-vanishing 
energy relaxation rate to the problem of dcphasing saturation observed 
at low T in the magnetoresistance of nanoscopic wires [4, 5, 2]. 

The infrared singular behavior of V(e), phenomenologically deduced 
above, means in particular that the interaction has no essential momen- 
tum dependence and should be of a local origin, in contrast to an inter- 
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action mediated by a diffusive or dispersive collective mode. Therefore, 
it has been proposed [6] that the anomalous energy relaxation may be in- 
duced by Kondo impurities, in particular by (nearly) degenerate atomic 
two-level systems [6] or other dynamical defects in the nanowire which 
can generate a two-channel Kondo (2CK) effect [7, 8]. The 2CK effect is 
known to have ideally a non-vanishing zeropoint entropy S(0) = A;filn\/2 
[9, 10] and, hence, can cause dephasing at low energies [11]. In fact, the 
electron interaction vertex mediated by a Kondo impurity does show a 
\je energy dependence. This was noted for the 2CK effect in the strong 
coupling region [6] (Fig. la)) and independently for the weak coupling 
(i.e. high energy) regime of the magnetic, single-channel Kondo (1CK) 
effect [12], see also [13]. Although this analysis takes elevated single- 
particle energies into account, it assumes the Kondo defect to be in 
thermodynamic equilibrium. Since, however, the effective interaction is 
generated dynamically due to a Fermi edge singularity, it will be substan- 
tially modified in the non-equilibrium situation of the energy relaxation 
measurements considered here. 

Although the Kondo effect in stationary non-equilibrium has been 
treated by various numerical methods before, no complete understand- 
ing has been reached up to now. In particular, the scaling properties and 
to what extent the Kondo effect at finite bias is a strong coupling problem 
have remained unclear. In this article we, therefore, summarize a pre- 
dominantly analytical analysis of the 1CK and the 2CK non-equilibrium 
Kondo effect [14]. It is shown that in nanowires with either 1CK or 2CK 
defects the quasiparticle distribution function exhibits scaling in terms 
of the applied bias, eU, when eU exceeds an intrinsic low-energy scale 
eU Q which, up to logarithmic corrections, is proportional to the equilib- 
rium Kondo temperature The numerical solution of the quantum 
Boltzmann equation gives quantitative agreement with the experimen- 
tal results discussed above, where in the scaling regime the impurity 
concentration is the only adjustable parameter. 

2. A SINGLE KONDO IMPURITY IN A 
NON-EQUILIBRIUM WIRE 

The 1CK as well as the 2CK system can be described by a generalized 
SU(N)xSU(M) Anderson impurity model in the Kondo limit, i.e. in 
terms of a quantum degree of freedom in a low- lying level e^, coupled 
to the continuum of conduction electrons via a hybridization matrix 
element V, where N = 2 denotes the (pseudo)spin degeneracy and M = 
1,2 the number of gegenerate, conserved conduction electron channels. 
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Figure 1 a) Leading contributions to the effective electron-electron vertex induced 
by a Kondo impurity, b) Diagrammatic representation of the NCA auxiliary particle 
selfenergies, Eqs. (1.5)-(1.7). Solid, dashed, and wiggly lines correspond to conduction 
electron, pseudofermion, and slave boson propagators, respectively (see text) . 



In the auxiliary particle representation [15] the model hamiltonian reads 

H = J2 e P c lma c pma + e d J2 fUa + V J2 Ul h ™ C P ma + h.C.) , (1.3) 

p,m,a a p,m,a 

where cj, mo . is the creation operator for an electron with momentum p. 

The auxiliary fermion and boson operators, fl, 6^, create the local defect 
in its quantum state a or in the unoccupied state, respectively. The 
operator constraint Q = J2a fafa+J2m "nfim = 1 ensures that the defect 
is in exactly one of its quantum states at any instant of time. In the 1CK 
case of a magnetic Anderson impurity a denotes spin, and m = 1 has no 
relevance. For a 2CK defect, a is identified with a pseudospin, e.g. the 
parity of the local defect wave function, and m = 1, 2 is the conduction 
electron spin acting as the conserved channel degree of freedom. Note 
that in this case the bosonic operator 6^ transforms according to the 
conjugate representation of SU(M) in order for the channel degree of 
freedom to be conserved. The equilibrium Kondo temperature of the 
model is (to leading exponential order) ~ Wexp(—1/NN J), with 
J = V 2 /|ed| the effective spin exchange coupling and W the half band 
width. 

The auxiliary particle representation allows for the application of stan- 
dard Green's function methods and, in particular, is readily generalized 
[16] to non-equilibrium by means of the Keldysh technique. It is known 
[17] to describe the universal spectral properties of the multi-channel 
Kondo effect correctly both in the infrared and in the high energy regime 
already within the non-crossing approximation (NCA), which is a self- 
consistent, conserving approximation to leading order in the effective 
hybridization T = V 2 ttN . Thus, in the multichannel case the NCA 
constitutes a correct infinite resummation of the logarithmic terms of 
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perturbation theory. In the 1CK case at finite bias the logarithmic se- 
ries is the same as for the 2CK effect in equilibrium because of the 
presence of two Fermi edges in the former case [14]. Thus, the NCA 
gives a correct description both of the 1CK and of the 2CK effects in a 
strong non-equilibrium situation (eU > and will be used for the 

explicit calculations below. 

In the remainder of this section we will consider a single Kondo defect 
at position x in a diffusive wire of length L with a bias voltage U applied 
between its ends. The generalization to the experimental situation of a 
dilute ensemble of Kondo impurities will be treated in the following 
section. The derivations apply both to the single and to the multi- 
channel case. For purely elastic scattering in the wire the electronic 
distribution function at the position of the impurity is 

f x (E, U) = jf°(E + eU) + (l - j ) f°(E) . (1.4) 

Denoting the momentum integrated greater (+) and lesser (— ) Keldysh 
Green's functions for the conduction electrons at position x by G+ X (E) = 
E p G+ x (E,p) = 2mf x (E,U) N (E), G~ X (E) = J2 p G~ x (E,p) = 
—2m(l — f x (E,U)) N (E), and similarly the "greater", "lesser" and re- 
tarded auxiliary particle propagators by Gj' b ' t (uj), the non-equilibrium 
NCA equations then read (Fig. 1 b)) 

w^w = -W^ G - ( - e)G ^ +e) (L5) 

qtM - — r de 

This set of non-linear equations is closed by the Kramers-Kroenig re- 
lations, G r j b (uj) = j de / (2-7ri) Gj b (uj) / — e + iO), which follow from 
analyticity and the fact that the auxiliary particle Green's functions 
have only forward in time propagating parts. The conduction electron 
t-matrix due to the Kondo impurity is 

It is crucial to determine the low-energy scales of the problem in 
non-equilibrium. These may be determined from perturbation theory. 
At bias eU the Kondo temperature, definied as the breakdown scale 
of perturbation theory, is suppressed, e.g. in the middle of the wire 
(x = 1/2), as 



= + N rj^-Gl x (e)Gf(u; + e). (1.6) 
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In addition, a Korringa-like spin relaxation rate l/r s oc eU appears, 
because in non-equilibrium there is finite phase space available for scat- 
tering even at T = 0, 

= }-NM (N J) 2 eU<eU . (1.9) 
T s (eU ) 47r 

Because of this inelastic relaxation rate, the logarithmic singularities of 
perturbation theory (T = 0) are shifted by l/r s to the complex fre- 
quency plane, and the Kondo scale disappears eventually for large bias 
eU and is replaced by the inelastic rate l/r s : The low-temperature scale 
is T (eU) = max[Tft-(eC7), l/T s (eU)]. The crossover from the Kondo to 
the Korringa scale occurs at a bias 

* - -hliM ^ - ^ ^ K|>) ■ <™> 

This means that in the large bias regime, eU > eU , the low-energy scale 
of the problem is proportional to the external bias itself. 

We now turn to the scaling analysis of the non-equilibrium Kondo 
equations (1.5)-(1.7). It is well known [18, 19], that at T = in equilib- 
rium the exact auxiliary particle propagators exhibit infrared powerlaw 
divergencies Gj b (co) oc 0(o>) u>~ a i' b , where the exponents q/^ are due to 
an orthogonality catastrophy between the occupied and the unoccupied 
impurity state. In the 2CK case the exact exponents are reproduced by 
NCA [17] and depend in a characteristic way on the spin degeneracy N 
and the channel number M, a f = M/(N + M), a b = N/(N + M) with 
ctf + ctb = 1. In the large bias regime (eU > eU ), the NCA equations are 
valid both in the 1CK and in the 2CK mentioned above. The 

behavior around the Fermi edges may be determined by defining the 
complex frequency variables z = oj — i/2T s (eU) etc. and rewriting the 
NCA equations for z —>■ and z — > eU as coupled differential equations 
in analogy to the equilibrium case [18]. We obtain damped powerlaw be- 
havior of the auxiliary fermion propagator for energies |cj| ~ l/r s and of 
the auxiliary boson propagator for energies around the two Fermi edges 
\lo\ ~ 1/t s , \uj — eU\ ~ 1/t s where the damping constant is l/r s . The 
exponents are modified in the large bias regime to a'j = 2M/(N + 2M), 
a' b = N/(N + 2M), where the additional factor of 2 originates from the 
fact that there are two separated Fermi edges. The form of a'p a' b is 
reminiscent of an effective doubling of the channel number M induced 
by the two Fermi edges. This behavior is confirmed by the numerical 
evaluation of the problem and is consistent with a recent perturbative 
renormalization group analysis [20] of the problem. The fact that power- 
law behavior (although damped) instead of logarithmic behavior persists 
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even at large bias may be seen as a remnant of the strong coupling be- 
havior of the Kondo problem in non-equilibrium. It is crucial for scaling 
in terms of ell to occur: Putting the powerlaw behavior of the auxiliary 
propagators back into Eqs. (1.5)-(1.7), using a'f + a' b = 1 and the im- 
portant fact that the low-energy scale T Q and the damping rate of the 
powerlaw behavior themselves scale with the external bias, T Q l/r s oc ell, 
it follows immediately that these equations depend only on dimension- 
less energies to/eU, E/eU etc., i.e. they are scale invariant in the regime 
eU > eU Q = 2NJ^j k B T^ m(-^y) (Eq. 1.10). It is seen that in 

experiments one may extract the equilibrium Kondo temperature 
from the lower breakdown scale of the scaling behavior. 

3. SOLVING THE KINETIC EQUATION 

The remaining task is to calculate the electronic non-equilibrium dis- 
tribution function f x (E,U) explicitly at an arbitrary position x along 
the wire in the presence of a dilute but finite density of 1CK or 2CK 

defects C'ifYip . 

The Boltzmann quantum kinetic equation for the momentum p depen- 
dent "greater" conduction electron Green's function G^ x ^(E,p), which 
contains information about the distribution function, reads [21] 

Ft + Vx " m + 6jg(f) ' V "HW^ = C ({ G W^}>^)> 

(1.11) 

where C is the collision integral induced by inelastic scattering off Kondo 
defects and E(x) are the external electric fields, including those gener- 
ated by static, random impurities in the wire. Exploiting the fact that 
the current in the disordered system is diffusive, 



1(E) = J2 L G+ x (E,p) = -DV xPx (E) = -DV x J2Gt x (E,p) 



p_ 

m 

(1.12) 

with D the diffusion coefficient, and summing in Eq. (1.11) over quasi- 
particle momenta p, the electric field term drops out as a surface term, 
and we obtain in the stationary case a diffusive kinetic equation [3] for 
the impurity averaged distribution function f x (E,U) = J2 p G^ x (E,p)/ 
2mN as function of the quasiparticle energy E and position x, 

-DV 2 J X (E, U) = C({J X (E, [/)}, c imp ) . (1.13) 

After momentum integration, the collision integral is given in terms 
of the impurity t-matrix (Eq. 1.7) to arbitrary order in the electron- 
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E/eU 

Figure 2 Scaled distribution functions f x (E,U) in nanowires. — Left panel: Two 
Cu nanowires (L = 1.5//m, D = 65cm 2 /s, x/L = 0.5 and L — 5.0^im, D = 45cm 2 /s, 
x/L = 0.5). Black curves: Experimental data for 0.05mV< U <0.3mV in steps of 
0.05mV [1]. The two curves showing deviations from scaling are at small voltages. 
Light curves: Theory for 2CK impurities in the scaling regime, eU > eU - Fitted 
2CK impurity concentration: Ci mp ~ 8 ■ 10~ 6 /(lattice unit cell). — Right panel: Two 
Au nanowires (wire 1 and wire 2 of Ref. [2]), measured at different positions x/L. 
Open circles: Experimental data; Solid lines: Theory for magnetic impurities in the 
scaling regime. Fitted magnetic impurity concentration: c imp « 2 • 10~ 4 /(lattice unit 
cell). 

impurity coupling T by 
C({/*(£,tf)},C; mp ) = |^^^ . (1.14) 

We have solved the set of non- linear integro-differential equations (1.5)- 
(1.7), (1.13), (1-14) numerically with the boundary conditions f x=0 (E, U) 
= f°(E), J X=L (E, U) = f°(E + eU) and that the electron system should 
be in local equilibrium with the dynamical impurity. The results are 
shown in Fig. 2. The theoretical f x (E/eU) curves show scaling behav- 
ior in the regime eU > ell (Eq. (1.10)) for bias voltages U varying 
within a factor of 4 to 9, depending on the high energy parameters of 
the model Eq. (1.3). There is excellent quantitative agreement between 
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theory and experiment for all samples. The experimental data for Cu 
and Au wires can be fitted equally well by the 1CK or by the 2CK 
theory, with very similar defect concentrations in the two cases. As an 
example, Fig. 2 shows the fits of the Cu data to the 2CK theory and of 
the Au data to the 1CK theory only. The agreement is consistent with 
the conjecture that the anomalous energy relaxation in those wires is 
mediated either by 1CK or by 2CK defects . 

4. RELATION TO DEPHASING 

For an ideal 2CK system in equilibrium, the non-vanishing quasipar- 
ticle scattering rate l/r(E) crosses over to a pure dephasing rate 1/t v , 
as the quasiparticle energy approaches the Fermi energy, E, T — > [11]. 
However, because of the finite level splitting A in real systems, one 
expects an upturn of at the lowest T ~ A//cb, with an interme- 
diate plateau for A/ks ~ T ~ T$ [11]. Such behavior can also be 
expected from magnetic (1CK) Kondo impurities in combination with 
a phonon contribution, because the spin flip rate reaches a maximum 
at T ~ [22]. This has been observed in the Kondo system AuFe 
[23] and most recently in clean Au wires presumably containing a small 
concentration of Fe impurities [2]. One might, therefore, conjecture that 
1CK or 2CK defects could be the origin of the dephasing time satura- 
tion observed in magnetotransport measurements of weak localization 
[4, 5, 2]. This assumption is indeed supported by several coincidences 
between the dephasing time measurements and the results on the non- 
equilibrium distribution function: (1) The dephasing time r v extracted 
from magnetotransport experiments [5, 4] is strongly material, sample, 
and preparation dependent. This suggests a non-universal dephasing 
mechanism, like dynamical defects, which is not inherent to the electron 
gas. (2) The dephasing time in Au wires is generically shorter than in Cu 
wires [5]. This is consistent with the fact that the estimates for the dy- 
namical impurity concentration Cj mp , obtained from the fit of the present 
theory to the experimental distribution functions, is substantially higher 
in Au than in Cu wires (Fig. 2). (3) In Ag wires one observes neither 
dephasing saturation nor E/eU scaling of the distribution function [24]. 
This is consistent with the assumption that dephasing saturation and 
anomalous energy relaxation in the nanoscopic wires have the same ori- 
gin and that there are no Kondo defects present in the Ag samples. 
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5. CONCLUDING REMARKS 

We have analysed the single- as well as the two-channel Kondo ef- 
fect in a stationary non-equilibrium situation. It was found that a 
Korringa-like, inelastic spin relaxation rate appears which at large bias 
sets the low-energy scale of the problem. Nevertheless, a remnant of the 
strong coupling Kondo fixed point persist even at large bias, which man- 
ifests itself by damped powerlaw behavior of the local spectral density. 
The latter leads to scaling of the non-equilibrium distribution function 
f x (E,U) = f x {E/eU) as an experimentally observable signature. The 
present theory yields quantitative agreement with the experimental re- 
sults for all samples measured, the density of Kondo defects in the wire 
being the only adjustable parameter in the scaling regime. This strongly 
suggests that the anomalous energy relaxation might be caused by either 
1CK or 2CK defects. However, it was shown that the scaling property 
of the distribution function does not distinguish between 1CK or 2CK 
impurities. Further experimental tests, like application of a magnetic 
field, will be required for that purpose. 
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